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Abstract 

Stationary four-dimensional BPS solutions to gravity coupled bosonic theories admit- 
ting a three-dimensional sigma-model representation on coset spaces are interpreted as null 
geodesies of the target manifold equipped with a certain number of harmonic maps. For 
asymptotically flat (or Taub-NUT) space-times such geodesies can be directly parametrized 
in terms of charges saturating the Bogomol'nyi-Gibbons-Hull bound, and classified ac- 
cording to the structure of related coset matrices. We investigate in detail the "dilaton- 
axion gravity" with one vector field, and show that in the space of BPS solutions an 
5*0(1,2) x SO(2) classical symmetry is acting. Within the present formalism the most 
general multicenter (IWP/Taub-NUT dyon) solutions are derived in a simple way. We also 
discover a large new class of asymptotically flat solutions for which the dilaton and axion 
charges are constrained only by the BPS bound. The string metrics for these solutions are 
generically regular. Both the IWP class and the new class contain massless solutions. 
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1 Introduction 



A typical feature of many non-linear field theories is the existence of certain positivity bounds 
for the energy similar to the Bogomornyi-Prasad-Sommerfield ( BPS) bound [[l| in the Yang- 
Mills-Higgs models admitting monopoles. Classical solutions saturating these bounds (BPS 
solutions) possess topological charges and minimize the energy for given values of these charges. 
Consequently they are stable, and usually associated with solitons. Physically the BPS satura- 
tion condition means that forces between them are balanced, so generically the corresponding 
multisoliton solutions also exist. In supersymmetric embeddings solitons possess unbroken su- 
persymmetries which are manifest in the existence of Killing spinors constant at spatial 
infinity. The vanishing of the corresponding spinorial variations on soliton backgrounds gives 
rise to linear Bogomol'nyi type equations which facilitate the effective construction of these 
solitons. Topological charges enter as central charges into the supersymmetry algebra and 
hence the knowledge of the soliton spectrum is important for a proper understanding of the 
dynamics of the theory. 

Similar phenomena are also well known in theories including gravity. The Einstein the- 
ory can be embedded in simple supergravity, the corresponding positivity argument providing 
a concise proof of the positive mass conjecture |J. The importance of BPS bounds in ex- 
tended supergravities was advocated by Gibbons 0. In N = 2 supergravity (containing in 
the bosonic sector the graviton and an Abelian vector field) the extreme Reissner-Nordstrom 
black hole saturates the BPS bound and has an unbroken supersymmetry, while the general 
non-singular solution of the corresponding spinorial equation is the Majumdar-Papapetrou 
multicenter solution to the Einstein-Maxwell (EM) equations ||. Since quantum corrections 
to supersymmetric backgrounds are controllable, the BPS solutions are an important tool to 
investigate the underlying quantum theory non-perturbatively [||. 

The N = 2 theory, however, has anomalies, so the non-renormalization argument is not 
fully applicable. The situation is better in N = 4 four-dimensional supergravity, containing in 
the bosonic sector the graviton, two scalar fields (dilaton and axion) and six U(l) vectors. The 
BPS bound for static asymptotically flat solutions to this theory was discussed by Gibbons 
Q. Truncated models of this kind containing one (dilaton-axion gravity), or two vector fields 
were extensively studied recently by Kallosh and collaborators M 0. In JiO| ] some stationary 
BPS solutions to dilaton-axion gravity were given and the BPS condition was generalized 
to include the NUT charge. Recently such solutions have been investigated in the context of 
heterotic string theory, which gives rise in four dimensions to a more general bosonic lagrangian 
including 28 vector fields, the dilaton and the axion and 132 scalar moduli fields forming a 
coset SO{6,22)/S(0{6) x 0(22)). The BPS bound in this theory was discussed in some 
static spherically symmetric solutions saturating this bound were studied in p2[ . Of particular 
interest are massless BPS states, which can substantially affect the low-energy dynamics fl3| . 

In all listed theories there are internal symmetries relating configurations depending ef- 
fectively on three space-time coordinates (e.g. static or stationary), which can be used to 
generate new explicit solutions. This was one of the main technical tools in most of the above 
papers. Another widely used approach was dimensional reduction ^from ten-dimensional so- 
lutions known to be exact string backgrounds, such as stringy plane waves fllO|| . Both methods 
generically give physically interesting, but nevertheless restricted classes of solutions. A more 
general approach consists in the direct use of the first order equations for the supercovariantly 
constant spinors. Along these lines Tod was able to find all stationary BPS solutions to EM 
theory, dilaton gravity, and a large (although restricted) class of solutions to A = 4 supergrav- 
ity [14]. However, the complete enumeration of all stationary asymptotically flat solutions to 
dilaton-axion gravity (and larger theories including moduli fields) seems not to have been given 
so far. Here we make some progress in this direction using an alternative (purely bosonic) way 
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to classify the BPS solutions. 

The method utilizes the cr-model structure of the three-dimensional theories describing 
four-dimensional backgrounds admitting a non-null Killing vector. For pure gravity this is the 
well-known Ernst SL(2, R) / SO(2) u-model [15], generalized to Maison's SL(n — 2,R)/SO(n — 
2) er-model [16] for n-dimensional vacuum Kaluza-Klein (KK) theory, while for the EM theory 
this is the Neugebauer and Kramer SU(2,1)/S(U(2) x U(l)) model 0, @. More general 
classes of such models following from supergravities were studied by Breitenlohner, Maison 
and Gibbons IS]. The particular case of dilaton-axion gravity (with one vector field) was 
studied in f20[| , where the coset structure Sp(A, R)/U(2) of the corresponding target space 
was found. Typically, when the Killing vector, with respect to which a reduction to three 
dimensions is performed, is timelike, the resulting target space has an indefinite metric and 
contains null geodesic lines. [] It is well known that geodesies of the target space equipped 
with some harmonic functions on a three-space generate a solution to the <7-model equations 
fl7|] , |pl]] . It was observed in [g2| that null geodesies of the target space of stationary five- 
dimensional KK theory may be used to generate multisoliton solutions similar to the Israel- 
Wilson-Perjes solutions |23j of EM theory. The association of stationary multisoliton systems 
with null geodesies in target space is a general feature of a-models, due to the fact that in the 
null case the three-space is conformally flat. This property can be related to the BPS force 
balance condition. Therefore, the classification of stationary BPS solutions may be traced to 
the investigation of the null geodesic structure of the three-dimensional u-model target space. 

We discuss the general features of this approach in Sect. 2, and then consider in detail 
the stationary dilaton-axion gravity model in Sect. 3. It is shown that isotropic geodesies 
in the space Sp{A,R)/U{2) fall into two classes depending on whether the the correspond- 
ing generating matrix is degenerate or not. The space of classical BPS solutions exhibits a 
symmetry under ,50(1,2) x S0{2) transformations mapping one asymptotically Taub-NUT 
geodesic solution onto another (Sect. 4). These transformations include electric-magnetic du- 
ality, mass - NUT duality as well as two special boost transformations in the non-degenerate 
case. Solutions related to degenerate matrices are shown to reproduce all previously known 
BPS solutions including rotating dilaton-axion Taub-NUT dyon [^] and dilaton-axion IWP 
solutions [1C], generalized to include all possible charges explicitly and to allow for arbitrary 
directions of rotation axes. Although the Einstein metrics for these extreme black hole solu- 
tions are singular, the corresponding string metrics are generically regular in the non-rotating 
case. Massless states with non-zero NUT-s are found. Non-degenerate matrices correspond to 
formally quasi-periodic solutions (Sect. 6). These solutions have dilaton and axion charges 
constrained only by the BPS condition. Their Einstein metrics present weak naked singular- 
ities, while the associated string metrics are again generically regular. Here we present these 
entirely new solutions in some detail, in view of their amazing simplicity. In Appendix A we 
describe a useful decomposition of the symplectic algebra, and find some intrinsic connection 
with 1+2 Clifford algebras. In Appendix B the IWP solutions to the EM theory are rederived 
within the same approach. 



2 General formalism 

Consider a four-dimensional action 

S = J (-^- + L(A I ,cf> a )) V=9d 4 x, (2.1) 

^^Note that the canonical procedure gives rise in the case of dilaton-axion gravity to the coset Sp(4, R) /U(l, 1) 
Jlif , but the indefinite metric can also be presented in terms of symmetric matrices belonging to Sp(4, R)/U(2) 
which are somewhat simpler. In what follows we assume this representation of the non-compact target space of 
stationary dilaton-axion gravity. 
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where A denotes the set of U{1) (linear) vector fields, and 4> a is the set of scalar fields typically 
including the dilaton, Peccei-Quinn axion and moduli. In extended supergravities/superstring 
effective theories scalar fields usually form coset spaces, while vector fields transform under 
a representation of the corresponding global symmetry group. When the theory is reduced 
to three dimensions, vectors can be traded for pairs of scalar fields, which can sometimes be 
incorporated into a larger coset space. The corresponding three-dimensional theory then turns 
out to be a gravity coupled sigma-model on a symmetric space. A general discussion and 
a list of different group combinations ensuring such a property were given by Julia [26] and 



Breitenlohner, Maison and Gibbons [19] and discussed later in a number of papers. Assuming 
a certain familiarity with this formalism, we just briefly recall here the basic formulas with an 
emphasis on the so-called geodesic solutions. 

We will be interested in stationary solutions depending on three (space-like) coordinates, or, 
in other words, admitting an (asymptotically) timelike Killing vector field. Then the standard 
dimensional reduction ansatz 

9 " U = {-M -/-%!+ MoJ (2,2) 

leads to the corresponding three-dimensional theory. In the stationary case, U{\) vector fields 
may be expressed in terms of electric v 1 and magnetic u 1 potentials, while the rotation one-form 
lu = ojidx 1 can be dualized to give the twist potential Xi the pair (/, x) forming a coset space 
SL(2, R)/SO{2) |p7| . To make u 1 and x dynamical variables one has to introduce Lagrange 
multipliers ensuring the Bianchi identities. The resulting action will be that of the gravity 
coupled three-dimensional c-model 

Sa = ^J (n- G A BdiX A djX B h ij ) Vhd 3 x, (2.3) 

where 1Z is the Ricci scalar built out of hij, and the set of scalars X A = (/, x, v 1 , u 1 , 4> a ), 
combined into the potential space endowed with the metric Gab(X), acts as a source of three- 
dimensional Euclidean gravity. The equations of motion for the X A 

1 „ / „ „ B 



^di [^hh^G AB d 3 X a ) = (2.4) 

define a harmonic map from the three-space {x 1 } to the potential space {A'" 4 } (target space). 

Obviously, the equations of motion ( |2.4j ) are preserved by the isometry group G of the target 
manifold. Therefore this group is a symmetry group of the space of solutions, using which new 
solutions may be generated. In the class of theories under consideration this group is large 
enough to ensure the symmetric space property of the target space, in which case the latter 
may be regarded as a coset space G/H, with H a subgroup of G. This implies the existence of 
a matrix M (an appropriate coset representative) in terms of which the target space metric is 
given by 

dl 2 = Gab dX A dX B = -^Tr (dMdAf -1 ) . (2.5) 
Consequently, the three-dimensional action will read 

Sa = ^J (n+ ^Tr (VMVM -1 ) j Vhd 3 x, (2.6) 

where V stands for the three-dimensional covariant derivative, and a scalar product with 
respect to the metric hij is understood. The corresponding equations of motion take the form 
of the conservation equation for the matrix current 

V (A/^Vm) = 0, (2.7) 
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while the three-dimensional Einstein equations read 



Kij = --Tr (ViMV^M^ 1 ) . (2. 



As was noticed by Neugebauer and Kramer [ |17| , when one makes the special assumption 
that all target space coordinates X A depend on x % through only one scalar potential, i.e. 
X A = X A [a{x 1 )]^ it follows from the equation of motion Q2.4D that this potential can be chosen 
to be harmonic, 

Act = 0, A = V 2 , (2.9) 
Eq. Q2.4| ) reducing then to the geodesic equation on the target space 

d 2 X A . dX B dX c 

^ + r «^r^r = °- (2 - 10) 

Rewriting this in matrix terms, 

*(M-1™)=0, (2.11) 



da V da 

one may express the solution to the geodesic equation in the exponential form 

M = Ae Ba , (2.12) 

where A £ G/H and B 6 £ie(G) are constant matrices such that M 6 G/H. Usually we 
are interested in asymptotically flat (Taub-NUT) solutions, and it is assumed that a{oo) = 0. 
Then the asymptotic value of M is A, some fixed quantity depending on the choice of the 
coset representative. Asymptotically flat solutions of this kind are thus target space geodesies 
passing through the point A. 



The three-dimensional Einstein equations (2.8) read, in the case of geodesic solutions, 

TZij = ^Tt(B 2 )ViaVja. (2.13) 
From this expression it is clear that in the particular case 

Tr(B 2 ) = (2.14) 
the three-space is Ricci-flat. In three dimensions the Riemann tensor is then also zero, and 



consequently the three-space is flat. From the Eq. (|2.5| ) one can see that the condition ( 2.14 ) 



corresponds to null geodesies |22[ of the target space 

dl 2 = ^Tr{B 2 ) da 2 = 0. (2.15) 

For different particular theories it was observed that null geodesies describe extremal black 
holes satisfying a Bogomol'nyi type condition on the relevant charges. The flat nature of 
the three space supports the anticipation that the attractive/repulsive forces associated with 
different charges are mutually balanced. In cases when the bosonic action can be embedded 
into some supergravity, it was argued that this condition corresponds to the existence of Killing 
spinors P] jl4|] ensuring the unbroken supersymmetries of bosonic solutions (usually called BPS 
solutions). Although no general proof is known so far that any BPS solution can be presented 
as a null geodesic passing through A, at least the inverse statement seems to be rather plausible 



m 



Whereas the group G exhibits a symmetry of the full solution space, only a restricted set 
of G-transformations preserves the class of asymptotically flat geodesic solutions. To find the 
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corresponding subgroup (in particular, we are interested in symmetries of the BPS subspace), 
consider a similarity transformation for the matrix B 

B^B V = UBU- 1 , U EG, (2.16) 

and construct the corresponding 

Mu = Ae Bua . (2.17) 

Since A is not transformed, Mjj is not necessarily an element of G/H. However for some U 
it may occur that Mjj £ G/H indeed, and thus is a solution. Such U will be the symmetry 
transformation mapping one geodesic solution onto another. The condition Tr(By) = will be 
satisfied automatically, if ( 233 ) holds. Therefore we deal with symmetries of the BPS space as 



well. 

The various null geodesic solutions fall into equivalence classes, corresponding to the in- 



equivalent types of matrices B satisfying ( 2.14 ). There is only one matrix type in the EM case, 



leading to the IWP class (see Appendix B). Three distinct matrix types are associated with null 



geodesies of the five-dimensional KK theory [22]. The first type yields regular quasi-periodic 



solutions, the regular solutions of the second type are extreme black holes, while the third class 
contains the KK magnetic monopoles and multipoles |24|| . As we shall see in the next section, 
only two matrix types occur in dilaton-axion gravity with one vector field. 

The construction ( 2.1 2| ) may be generalized |22[] to the case of several harmonic functions 



a- 



Aa a = 0, (2.18) 

by showing that 

M = AexpQ2 B * a *) ( 2 - 19 ) 



solves the field equations fl2.7| ) provided that the commutators [B a ,Bb] commute with the B c : 

[[B a ,B b ],B c ] = 0. (2.20) 
To prove this, we first rewrite ( |2.19| ) as 

M = Aexp(-i^^[ J B 6 ,B c ]a 6 a c ) JJe 8 -'- (2.21) 

c>b b a 

where the first exponential may be treated as a c-number function. We then obtain for the 
matrix current 

a/- 1 vm = Yu B * Va * -jEE^- B ^ fa v<7 <= - > ( 2 - 22 ) 

a, c>b b 



which is conserved if the a a are harmonic. The three-dimensional Einstein equations ( |2.8p 
generalize to 

Rij = lEE Tr(B a B b ) V l a a V J a b , (2.23) 



4 

a b 



so that the three-space is Ricci flat if the matrices B a satisfy 

Tr{B a B b ) = . (2.24) 
So the number of independent harmonic functions on which may depend a BPS solution of 



the form ( 2.1S| ) is limited by the number of independent mutually orthogonal null vectors of 



the target space. In the class of theories we are dealing with, the target space is a locally 
Minkowskian space M Pj(? , where p counts the positive eigenvalues coming from the two gravita- 
tional potentials /, x an d the effective scalar fields, and q is the number of negative eigenvalues 
coming from electric and magnetic potentials. Then the number of independent null vectors is 
inf (p, q). The actual number of possible independent potentials, however, may be less because 



of the extra condition ( p. 20 ). 
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3 Dilaton-axion gravity 



In the following we turn to a specific theory which may be regarded as a (truncated) bosonic 
sector of the four-dimensional heterotic string effective action, or that of the N = 4, D = 4 
supergravity with only one non-zero vector field. The model describes the gravity-coupled 
system of two scalar fields (dilaton 4> and (pseudoscalar) axion k), and one Abelian vector field 



A,: 



S = lib / {~ R + 2d ^ d ^ + l eA % Kd ^ - e-^F^F^ - kF^F^X ^d 4 x, (3.1) 

where F^ = \E^ vXt F Xt , F = dA . 

Assuming stationarity, one can perform reduction to three dimensions using the metric 
ansatz (|2.2|). The vector field can be parametrized by 2 real functions: an electric potential v, 



F iQ = d iV /V2, (3.2) 

and a magnetic potential u, 

e -2<t> F ij + K pij = ffijkQ^i^ (33) 

The corresponding spatial vector potential Ai may be recovered from u and v by 

V2dA = e 2 ^f~ 1 *(du-Kdv)-dvAio. (3.4) 

To satisfy the mixed components of the four-dimensional Einstein equations one has to 
express u>k through the twist potential x as follows 

Ti = d iX + vd iU - udiv, t* = -fe ijk djUJ k /Vh. (3.5) 

Here and below 3-indices are raised and lowered using the three-space (euclidean signature) 
metric hij and its inverse h lJ . 

The resulting three-dimensional gravity coupled cr-model has a six-dimensional target 
space, X A = (/, X-, v , u , 4>i K )i A = 1, 6 endowed with the metric 

11 1 

dl 2 = —2 {df 2 + (d X + vdu - udv) 2 } - - [e'^dv 2 + e 2(t> (du - ndv) 2 ) + 2# 2 + -e^dn 2 . 
2j j 2 

(3-6) 

A representation similar to ( |3.6[) has been found for the stationary EM system by Neuge- 



bauer and Kramer [17], their formula can be recovered by setting (j) = k = 0. As was shown by 
Mazur [p^| , the EM target space is isomorphic to the symmetric space SU(2, 1)/S(U(2) x U(l)). 
When the dilaton and axion fields are introduced one obtains instead the six-dimensional sym- 
metric space Sp(4, R)/U(2) |2(J. A matrix representative of this coset can be chosen to be the 
symmetric symplectic matrix [^7]] 

M= ( QP- 1 P+QP^Q ) ' (3 ' 7) 
where P and Q are the real symmetric 2x2 matrices 



with w = u — kv. 
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A solution M depending on only one potential is a geodesic 

M = Ae Ba (3.9) 

passibg through a given point A. If we choose as usual the harmonic potential to vanish at 
spatial infinity, then A = M(oo). For the asymptotically flat and/or asymptotically Taub- 
NUT configurations we are interested in, the "potentials" X A should be normalized at spatial 
infinity by 

/(oo) = 1, x(°°) = v(oo) = u(oo) = 0(oo) = k(oo) = 0, . (3.10) 



A=( ? ° ) (3-11) 



leading to 

A = ( 

(7 3 . 

(here and below / is the 2x2 unit matrix, o\ = a x , 02 = icr y , 03 = a z , and a x , a y , a z are the 
standard Pauli matrices with a z diagonal). Now we state the conditions on the matrix B such 
that M £ Sp(A, R)/U{2). The symmetry of M leads to the pseudosymmetry condition 

B T A = AB. (3.12) 

M is also symplectic, M T JM = J, with 

'=(-°, 0) 



K=[ » J). (3.15) 



On account of ( 3.12 ) and of the symplecticity of A, this leads to the condition 

BK + KB = 0, (3.14) 

with 

K=( 

-a 3 

Finally, M is unimodular, leading to 

TtB = 0. (3.16) 

As shown in Appendix A, the exponential e Ba belongs to the coset Sp(4, R)/(SO(2) x 
SO(l, 2)), where the SO (2) component is generated by K, while the remaining non-compact 
group is generated by the E , 

(note that So = J), with the commutation relations 

[S , E 6 ] = 2e a6 c E c , [K, E a ] = . (3.18) 

Here and below we use <SO(l, 2) indices a,b,c = 0, 1, 2 which are lowered/raised using the 1+2 
Minkowski metric rj a b = diag(l, —1,-1) and its inverse rj ab = 7]^. 

The six remaining elements of sp(4, R) can be cast into two so(l, 2) (co)vectors 

" "/). (S -Ml. (»») 



and 



satisfying the commutation relations ( |A.14[) , (|A.15 ) with and E a . 
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Therefore, the most general B can be parametrized by two SO(l, 2) vectors a, (3 as follows 

B = a a r x a + a T 2 a = a ■ T 1 + (3 ■ T 2 . (3.21) 
To perform exponentiation one needs an anticommutator 

{r 1 ,, r 2 ,} = -2e ab c rf d r° d , (3.22) 

which lies outside the Lie algebra sp(4, R): 

M(? :)• (0 c ?)}■ 

Remarkably, this third triplet of matrices makes the full system SO(l, 2) covariant on the upper 
index as well, as discussed in Appendix A. The anticommutator 

{r c a , r d b } = 2 (VV / - e cd f e ab e r/ e ) (3.24) 

is valid for all c, d. This last relation shows that the three sets of gamma's form 1+2 Clifford 
algebras. 

Using these formulas one immediately obtains for the square of ( |3.21| ) 

£ 2 = - (a 2 + /3 2 ) / - 2 (a A /3) • r°, (3.25) 
so that the null geodesic condition ( |2.14[ ) reads 

Tr(£ 2 ) = -4(ct 2 + 1 ) = 0. (3.26) 
Assuming this condition to hold, we obtain successively 

B 3 = -2 (a' • T 1 + ■ T 2 ) , (3.27) 



where 

a' = /3A (a A/3), = a A (/3 A a), (3.28) 

and 

B 4 = 4(aA/3) 2 /. (3.29) 

The fact that the matrix B 4 is proportional to unity could also have been derived from the 
matrix identity 

V(B) = (3.30) 

where 

4 

V(X) = det(B-XI) = Y[(X-Xi) (3.31) 

is the characteristic polynomial. Because the matrix B is constrained by equations (|3.16| ) and 
( |3T26[ ), TrB = and Tr(B 2 ) = 0, leading to Tr(B 3 ) = as well from eq. (|]27]), only zero and 
fourth powers of A enter this characteristic polynomial. Consequently, the identity ( 3.30 ) reads 

B A + detB = 0. (3.32) 

Now we are in a position to classify the possible types of matrices B satisfying the null 
geodesic condition fl3~2^ ). It follows from Eqs. ( |3~29| ), ([T26|) and ( |3~32D that there are two 
essentially different types of B matrices: 
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i) Type 1: degenerate B, det B = 0. 

This is the case whether la: a and (3 are collinear, 

(3 = cot, < c < oo, with a 2 = (3.33) 

^from Eq. (3.26), or lb: a = 0, 1 = or vice versa (Eq. (3.33) goes over to this in 
the limit c — > oo or c — > 0). For the type 1, .B 2 = and hence the expansion of the 
exponential exp(Ba) contains only a linear term. 

ii) Type 2: non-degenerate B, det B / 0. 

type 2a : a and (3 are non- collinear, but both are null, i.e. a 2 = 0, /3 2 = 0; 
type 2b : neither of a and (3 is null, that is either 

a 2 = -f3 2 > 0; (3.34) 

or 

(3 2 = -a 2 > 0. (3.35) 

As we will see later, there is an internal symmetry of the space of geodesic solutions mixing 
subclasses a and b inside each type, but preserving the distinction between types 1 and 2. 
Now we look for BPS solutions depending on two harmonic functions o~, r 

M = Ae Ba+CT , (3.36) 

with C satisfying the same conditions as B, and thus admitting an expansion 

C = 7-T 1 + (5-r 2 , (3.37) 

and, in addition^, 

[B, C] = 0, Tt(BC) = 0. (3.38) 
These last equations impose the following conditions on the SO(l,2) vectors 7, 6: 

a A 7 + (3 AS = 0, 
aS-p-y = 0, (3.39) 
a • 7 + (3 • S = 0. 

These conditions lead to very strong restrictions on possible 7, S once a, (3 are given. 
Namely, for B of type 1, the matrix C should be also of type 1, and 7, 8 should be collinear to 
a, (3; it then follows that BC = CB = 0. Redefining a, r through suitable linear combinations 
one can then always choose 

B = (0,/3), C=(/3,0), with /3 2 = 0. (3.40) 

Since for the type 1 the exponents in ( |3.36 ) contain only linear terms it is clear that the most 
general solution can be presented in the equivalent form 

m = a{i + (3- (rV + rV)}. (3.41) 

If B is of type 2, then the only solution is C oc B (i.e. there is no extreme solution really 
depending on two potentials). Since the target space ( |3.6| ) is locally M^, it is also clear that 
for the present system no BPS solutions depending on more than two potentials are possible. 

2 The study of the Lie algebra sp(4,R) shows that the more general condition ( |2.19| ) on the commutator [B, C] 
reduces to [B,C] — 0. 
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4 Physical charges and S0(1, 2) x 50(2) symmetry of BPS space 



From the above analysis it follows that any stationary asymptotically flat BPS solution to 
dilaton-axion gravity can be parametrized by six real numbers subject to the null-geodesic 
condition (|3.26| ). These numbers can be conveniently arranged into two S0(1, 2) vectors a, f3. 
Here we want to relate them to physical charges which can be introduced via the standard 
asymptotic expansions of the relevant fields in the asymptotically Taub-NUT space-time (in- 
cluding the usual asymptotically flat one as a particular case). Here we introduce the charges 
in accordance with the paper [25] as follows (recall the gauge conditions ( 3. 1Q| ) ) 



r 

2N 



X 



r 

V2Q 



r 

u ~ (4.i) 

r 

r 

2A 

k ~ — , 

r 

It is also convenient to use a complex mass combining the Schwarzschild mass M and the NUT 
charge N: 

m = M + iN, (4.2) 
a complex axidilaton charge combining the dilaton D and axion A charges 

d = D + iA, (4.3) 

and a complex electromagnetic charge 

q = Q + iP, (4.4) 

joining the electric Q and magnetic P charges. 

To make the identification of the components of the vectors a, f3 with physical charges, we 
consider a solutio n (|2.12D depending on a single monopole potential a, with the matrix B of 
the generic form ( |3.21| ). We choose the vectors cx, (3 to have the dimension of a length (i.e. the 
dimension of the charges in flf.ip ). The dimension of the harmonic function a will then be an 
inverse length, and we choose the normalization 

a ~ - . (4.5) 

r 



Substituting the expansions (|4.1| ) into the matrix M in the form (|3.7] ) , and comparing with the 
corresponding linearization of (|3.9|), one finds the following correspondence 



a = (V2P, A- N, M + D) , 

(3 = (-V2Q, M-D, N + A) . (4.6) 



The null geodesic (BPS) condition |T26|) a 2 + 1 = then assumes its standard form 



Q 2 + P 2 = M 2 + N 2 + D 2 + A 2 , (4.7) 
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or, in complex notation, 

\q\ 2 = \m\ 2 + \d\ 2 . (4.8) 

For type 1 and type 2a some additional constraints are imposed on a, (3. Consider first 
the degenerate type 1. The collinearity condition (without assuming the BPS constraint ([4.7])) 
gives 

M 2 + N 2 = D 2 + A 2 , (4.9) 

Q(N-A) = P(M-D), (4.10) 

P(N + A) = -Q(M + D) (4.11) 

(only two of these equations are independent). Solving these equations with respect to D, A 
one gets 

M(P 2 - Q 2 ) - 2NQP 
Q 2 + P 2 ' 
N{Q 2 -P 2 )-2MQP 

A ~ Q2 + P 2 ' ^ AZ > 



d=-™>. (4.13) 

q 



or, in complex form 

Now, if the BPS condition ( |4.7|) is imposed, one gets from (^. 

Q 2 + p 2 = 2 (M 2 + iV 2 ) , (4.14) 



or qq = 2mm, so that the relation (|4.13j ) takes the "holomorphic" form 



q 2 

d=- — . (4.15) 

2m v ; 



This relation was found for a rotating Taub-NUT dyon in dilaton-axion gravity in |25j. Par- 
ticular cases of these dyon configurations are the purely electric type P = (a = 0), for which 
the extremal states have dilaton and axion charges 

D = -M, A = N, (4.16) 

and the purely magnetic type Q = (/3 = 0), for which one has 

D = M, A = -N. (4.17) 

Let us now consider the possibility of solutions with zero Schwarzschild mass. Within the 
type 1 such non-trivial configurations may exist if the NUT parameter is non-zero. In the 
generic dyon case (type la) the dilaton and axion charges are 

QP Q^F! 

N 2N v ' 

and the BPS condition reads Q 2 + P 2 = 2N 2 (magnetic solutions of this kind were given in 
|31j| ). In non-dyon cases 

L> = 0, A = ±N, (4.19) 
where the upper sign coresponds to electric and lower to magnetic solutions. 
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Now we discuss the relations between charges for solutions corresponding to non-degenerate 
B. In the null case 2a one has in addition to the BPS condition (4.7) a single constraint 



AID - NA 

Contrary to type 1, in the case 2a generically 



P 2 -Q' 



M 2 + N 2 ^ A 2 + D 2 



(4.20) 



(4.21) 



(if the equality holds we come back to the case la). Solving Q4.7| ) and ( 4.20 ) we obtain for the 
dilaton and axion charges in the case 2a 



D 
A 



where 



Comparing this with the "standard' 



M (P 2 - Q 2 ) - 2NQPg 

2(M 2 + N 2 ) ' 
N (Q 2 - P 2 ) - 2MQPg 

2 (M 2 + N 2 ) ' 

[Q 2 + P 2 - 2(M 2 + N 2 )] 2 
4Q2p2 • 

relations (4.12) one can see their difference unless p 



(4.22) 



(4.23) 



1 



in which case we come back to the type 1. Finally, in the non-null non-degenerate type 2b 



case the physical charges are constrained only by the BPS condition (4.7). 

Now let us discuss the symmetry transformations acting inside the class of stationary asymp- 
totically Taub-NUT geodesic solutions to dilaton-axion gravity. As was explained in Sect. 2, 
we look for similarity transformations for the matrix B which preserve the pseudosymme- 
try condition ( 3.12| ). It is shown in the Appendix A that if M € Sp(A, R)/U(2) and A is 



given by ( pUp , then AM = e Ba £ Sp(4,R)/(SO(l,2) x SO(2)). Therefore the subgroup 
H' = SO(l,2) x SO{2) is the desired symmetry group. In algebraic terms, to preserve the 
condition ( |3. 12[) , the generator X E sp(4,R) of the similarity transformation ( p. 16 ) should 
satisfy 

X T A = -AX. (4.24) 
In view of the decomposition of the sp(4, R) algebra described in Appendix A, it is clear that 



X E TC (cf. (All)). In fact this subgroup of the symplectic group has a larger meaning as 
that preserving asymptotic flatness (so it is also relevant to non-null geodesic solutions). It is 
worth giving a more detailed description of this symmetry with an emphasis on its action on 
BPS solutions. 

The group H' consists of two components. The first one (SO (2)) is generated by K, i.e 
(12.16|) U = e Ke/2 . This transformat ion mixes a and /3: 



m 



ex — > a cos e + f3 sin e, 
/3 — ► /3 cos e — a sin e. 



(4.25) 



In terms of charges this is an electric-magnetic duality rotation, i.e. the transformed quantities 
are 

q' = qe~ ie , m' = me^, d! = de~ ie . (4.26) 

It is worth noting that when Q and P get mixed by this transformation, at the same time the 
mass and the NUT charge are also mixed, i.e. one is not allowed to make a dyon from a purely 
electric (purely magnetic) asymptotically flat (NUT-less) configuration without introducing a 



NUT charge. Clearly, the axidilaton charge follows the transformation rule implied by ( 4.15 ). 
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Comparing this with the general case of the action of the symplectic group presented as the set 



of scale, gauge, S-duality, Harrison and Ehlers transformations [27|, one can observe that the 
^-transformation is a certain combination of Ehlers, gravitational gauge and two S-duality 
transformations. 

Now consider the SO(\, 2) part of H' . It is generated by the S a accor ding to ( [Up . In 
view of ( A.15| ), these are the transformations acting independently (and similarly) on a and 



s e/2. 



(3. Pure 'spatial' rotations are given by U = e 

a -» a , 

a 1 -> a 1 cos + a 2 sin 0, (4.27) 
a 2 — ► a 2 cos — a 1 sin 
(similarly for /3). In terms of charges we find 

q' = q, 

m! = me"' 9 , (4.28) 
d! = de ie . 

This transformation can be used to compensate the mass - NUT charge mixing when the K— 
transformation is applied, if strict asymptotic flatness is desired to be preserved. On the other 
hand, this transformation may be used to generate a NUT charge for purely electric and purely 
magnetic configurations. It is also a (different) combination of Ehlers, gravitational gauge and 
two of S-duality transformations. 

The two other SO(l,2) transformations are boost-like, they have different actions on null 
and non-null a and (3. Consider first the non-null case. The boost in the plane — 1 is 
generated by U = e S2 ^ 2 : 



a 



o 



a° cosh £ — a 1 sinh £, 



a 



1 -> a 1 cosh £ - a sinh £, (4.29) 



2 2 

a — > a . 

Together with the corresponding transformation of j3 this results in 

Tfi — d 

q = gcosh^H ^^sinh£, 

v2 

m' + d' = m + d, (4.30) 
m! — d' = (m — d) cosh £ + V2q sinh £. 

This is a combination of electric gauge and Harrison transformations which change all six 
charges in a non-trivial way. Similarly, the — 2 boost U = e El7? / 2 gives 

q' = qcoshrj + i — ^^sinhr/, 
v2 

m'-d' = m -d, (4.31) 
m' + d' = (m + d) cosh r\ — iV2q sinh rj. 

This corresponds to a combination of magnetic gauge and Harrison transformations. 

If a and (3 are null, the boosts reduce to rotations and rescaling. Hence the action of the 
full SO(l, 2) component of H' is in this case two-parametric. If in some frame 

a = a°(l,n), (4.32) 
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where n is a unit two-dimensional vector, in any other frame one has 

ol = a'°(l,n'), (4.33) 

where n' is a rotated version of n. Note that (3 should be transformed using the same param- 
eters. Taking also into account the additional degree of freedom due to the action of K, we 
see that for a given potential a the general type 1 solution depends on three parameters (we 
do not count here for the asymptotic values of the dilaton and axion fields which can easily be 
induced once the solutions in the present gauge are given). 
Generally both a 2 and /3 2 as well as the scalar product 

a - (3 = -2{QP + AM + ND) (4.34) 

are invariant under the SO(l, 2) transformations. These quantities, however, are not preserved 
by the electric-magnetic duality transformation ( [4.25| ), so if one wishes to identify solutions 
related by symmetries, formally there is no distinction between subcases a and b in the clas- 
sification of the previous section. However the distinction still may be useful because the 
corresponding solutions are physically different. 
In the general case a rescaling 

ex — ► ken, (3 -> k(3, (4.35) 

is also an invariance transformation although trivial, equivalent to a redefinition of harmonic 
functions, physically it corresponds to an equal rescaling of all six charg change of 

length scale. 



5 Type 1 solutions 



All solutions belonging to the degenerate type 1 can be presented in the form ( 3.41| ), generally 
they depend on two real potentials a, r or one complex potential 

C = a-ir. (5.1) 

Contrary to the previous section, here we choose this potential to be dimensionless, so that 
the SO(l, 2) null vector (3 in ( 3.41 ) will also be dimensionless. Normalizing this vector so that 
(3° = 1, one can write 

(3 = (1, cos a, sin a). (5-2) 

Keeping in mind that from any solution one can generate by H' transformations a two- 
parameter family of solutions, as described in the previous section, we could start from a 
given value of a (say a = 0) in ([T^) . However we find it instructive to leave here the choice of 
the gauge a open. The corresponding solution ( 3.41] ) may be written in complex form 



M = ARe {l + 13 ■ (r 2 + iT 1 ) c} 



(5.3) 



Identifying this with the initial coset representation (3.7) we obtain the sigma-model vari- 
ables in terms of the harmonic potentials a and r: 



z = k + ie 



r 1 

X 

v + iu 
2d> 



l+<7 cos a + t sin a, 
/(rcosa — fjsina), 

-K, 

1 + C cos a 
1 + i^sina 



(5.4) 
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Comparing this with the corresponding form of the BPS solutions to the EM system (see 
Appendix B) one can observe that the scale factor f~ l now depends linearly on the harmonic 
functions, contrasting with a quadratic dependence in the EM case. As we shall see, this has 
a dramatic effect on the regularity property of the solutions. 

Let us consider some particular families of BPS solutions corresponding to different choices 
of the complex harmonic function £. 



Extremal rotating Taub NUT dyon 



We take the harmonic function £(r) to be 



c 



V2 



q 



\r + ia\ 



and choose spheroidal coordinates (r, 9, ip) such that [23| 

(r A a) 2 = a 2 (r 2 + a 2 ) sin 2 9 , r ■ a = ar cos I 



(5.5) 



(5.6) 



(a = |a|), and ip is the azimuthal angle around a. Identifying the physical charges by comparing 
Q4.6]) and (|5.2[), and taking into account the relation 



V2™ 



which results from ( 4.13| ), we obtain the solution 



ds 2 



A — a sin 2 



(dt + tvdip) 2 — X 



dr 



A sin 2 ■ 



x + d ° 2 + A 
A A — a z sin i 



■dip 2 



where 



A = r 2 + a 2 , 
S = r 2 + a 2 cos 2 + 2(iVa cos 9 + Mr), 
2 { (r 2 + a 2 )N cos 6> + aMr sin 2 9} 



The electric and magnetic potentials are 

V2 



r 2 + a 2 cos 2 
V2 



^ (Qr + Pa cos #) , n = (Pr — Qa cos 9) , 



and the complex axidilaton is 



r + m — d + ia cos i 
r + m + d + ia cos I 



The corresponding dilaton factor reads 



= 20 



1 

E 



r + 2M + ia cos 9 



Qq 



in 



(5.7) 



(5.8) 



(5.9) 



(5.10) 



(5.11) 



(5.12) 



This solution corresponds to the BPS limit of the general rotating Taub-NUT dyon solution 
found in |2f| (the radial coordinate used there differs from the present one by a constant shift 
r — > r — 2M). For a/0 the solution has a naked singularity, as any rotating BPS solution to 
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this theory in four dimensions [28]. For a = the metric has a simple form (previously given 
by Kallosh et al @) 



ds 2 = (1 + 2M/r)~ 1 ( dt + 2N cos 9 dip) 2 - (1 + 2M/r) dr 2 -r(r + 2M)(d8 2 + sin 2 6 dip 2 ) . (5.13) 
This extreme black hole metric turns out to be singular, the area of the would-be horizon at 



being zero. It is pointed out in pi that the monopole string metric ds z 



^ds 2 



IS 



regular for certain parameter values. Actually we have found that it is always regular, except 
in the lower-dimensional domain PM — QN = 0. The string metric corresponding to ( |5.13 ) is 



ds 2 



^ds 2 



(dt + 2N cos9 dip) 2 



(r + 2M)' 



with 



2/~>2 



r = (r + jpy + yq 



PM - QN 



(5.14) 



(5.15) 



Q 2 + P 2 " 

For Q 0, r is positive definite, and the "throat" r = is at infinite geodesic distance. In the 
purely magnetic case Q = 0, T = (r + 2M) 2 leading to the same conclusion, except if M = 
(then also D = 0, A = — N, P 2 = 2N 2 ) where the string metric is massless Taub-NUT. In the 
special case PM — QN = (7 = 0) with M/fl, the string metric near r = 



ds 2 ~ r 2 (dt + 2N cos 6 dip) 2 



dr' 



dW 



(5.16) 



has a horizon of zero area at the origin r = of Euclidean space. This "black-point" singularity 
1 29] is mild, all non-radial geodesies being deflected away from r = 0. 

The a / solution depends on the four parameters a, q (complex), and a. The first three 
parameters correspond to the three degrees of freedom of type 1 solutions under the action of 
the group 50(1,2) x 50(2), while the parameter a may be associated with the scale of the 
dipole potential r. However we should keep in mind that our monopole-dipole solution has 
been obtained from the choice ( 3.40 ) of the matrices B, C in (|3.36j ). A more general choice 
solving the constraints ( 3.3S| ) is, up to 50(1,2) x 50(2) transformations and rescalings of a 
and t, 

5 = (0,/3), C = {(3 cosip,(3 sinV) . (5.17) 
The resulting linear form of the matrix M 

M = A{l + /3 • (r x r cos^ + T 2 (a + r sin V))} (5.18) 

with the same potentials a and r as in Q5.5|) leads to the general monopole-dipole BPS solution 



depending on five parameters (the same solution could alternatively be obtained from ( 3.41 ) 
by transforming the potentials a — > a + r sin ip , r — > r cos ip). 



Multicenter solutions 



Let us first discuss the equilibrium conditions between any pair i, j of the multicenter system 
in which each center in endowed with mass, NUT charge, electric and magnetic charges, and 
dilaton and axion charges obeying to ( |4.7| ), ( ]4.9[ ) and ( 4.13j ). The condition that the various 
attractive (tensor and scalar) and repulsive (vector) forces balance is Eq.( 2,24| ), which translates 
into 

QiQj + P % Pj = MiMj + N t Nj + DiDj + AiAj (5.19) 
or, in terms of complex charges, 



Re {q%qj — rrnrrij — didj) = 0. 



(5.20) 
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Choosing here the general representation of the matrix Bi in terms of three parameters //j, «j, 

cxi = 2/Xj Ui cos ipi Pi = 2/ij Ui sin ^ , (5.21) 
where ui = (1, cos ctj, sin aj), Eq.( 5.19| ) reads 

on ■ cxj + f3 { ■ Pj = ifiifij cos(V>j - ipj) Ui -Uj = . (5.22) 

This condition has two possible solutions. The first solution is Ui ■ Uj = which, for null 
vectors itj, uy, means ay = a^, so that the ay, (3j are collinear to the tXj, /3j-; then the whole 
system of conditions ( 3.3S| ) is satisfied, so that BiBj = BjBi = 0, and 

M = A [1 + + Sj-oy] (5.23) 

is a two-center BPS solution, which can straightforwardly be generalized to any number of 
centers. The condition ( |5.22| ) is also solved by cos(V>i — ipj) = 0, or 

QiQj + PiPj = , (5.24) 

with Ui ■ Uj ^ 0; the other conditions ( 3.39| ) are however not satisfied, so that the equilibrium 
condition ( 5.24 ), which can hold only for a two-center system, does not allow for a geodesic 
BPS solution. 

Now we consider an equilibrium configuration of n centers r = rj endowed with complex 
masses rrij = Mj + iNj and charges qj = Qj + iPj (as well as induced axidilaton charges 
dj = —qj/2nij) subject to the condition that the angle parameter ay entering the matrix Bj, 
or (by virtue of ( |5.7|) ) the complex ratio m,j/qj, is independent of j. We also assume that 
each of these centers carries an arbitrarily oriented dipole moment Oj. Then, to write down 
the multi-rotating Taub-NUT dyon solution one merely has to change the complex harmonic 
function ( |5.5| ) to 

n 

C = -^E^' R) = {r-r 3 +ia 3 ) 2 . (5.25) 

This solution generalizes the axion-dilaton IWP solution presented by Kallosh et al [ 10 1 to 
the dyon case, presence of NUT's and arbitrary directions of rotation vectors for each center. 
Clearly, for aj / we deal with naked singularities. 
In the non-rotating case, the metric has the form 

ds 2 = L + p^)j U+p 2N 3 cosO, &p\ -(l+p^j dr\ (5.26) 

where dr 2 is the flat space line element, and 8j, ipj are polar angles relative to an arbitrarily 
oriented orthogonal frame centered at r = rj. An essential difference with the Majumdar- 
Papapetrou (MP) solution to the EM system (see Appendix B) is that the metric function 
/ now appears as an inverse linear function of Coulomb potentials, while the MP solution 
(endowed with NUT-s) 

ds 2 = \U\- 2 (dt + u- dr) 2 - \U\ 2 dr 2 , (5.27) 

with 

U = 1 + V j „ S , V Au = lm[2VU + (UVU -UVU)}, (5.28) 

has quadratic factors. This is the reason why the surfaces Rj = , which for the MP solutions 
are regular horizons [30], correspond in the dilatonic case to space-time singularities (the radius 
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of two-spheres shrinks to zero) . It was shown recently || that adding a second vector field one 
can get dilatonic solutions sharing the above MP property too. The other difference between 
flUD and ( ET27D is of course the non-linear dependence ( |5.28| ), in the EM case, of the rotation 
one-form uj on the iVj and Mi (except if all the iV, or all the Mj vanish); this is a consequence 
of Eq.dpD for the non-commuting matrices B and C of Appendix b. 



Massless states 



BPS solutions with a vanishing Schwarzschild mass are of particular interest in the under- 



lying quantum theory since they can substantially influence the low-energy dynamics [13|. As 
was noticed in the previous section, the type 1 solutions may be massless but at the expense 
of a non-zero NUT charge (so that the complex mass m remains finite), the BPS condition 
now reading Q 2 + P 2 = 2N 2 . Although they are not asymptotically flat in the usual sense, it 
is worthwhile giving these solutions in explicit form. All non-rotating massless solutions turn 
out to possess a locally flat metric 



ds 2 = (dt + 2N cos 9 dip) 2 - dr 2 



\d9 



2 + sin 2 9dtp 2 ) 



(5.29) 



which is massless Taub-NUT. The dilaton factor is asymmetric with respect to magnetic and 
electric components: 

9rh 2QP 2Q 2 

so for purely magnetic solutions = 0, while for purely electric ones the dilaton field is non- 
zero, but short-range. In the general dyon case the dilaton charge is non-zero. The axion 
field 

2(Ar-ND) 

K ~ ( r + D )2 + (jv + A) 2 V- AL ) 

is non-zero in all cases, the values of the dilaton and axion charges being given by (4.18) or 
( 4.19 ). In non-dyon cases the axion field is short-range. 

Massless multicenter solutions may be constructed in a similar way. In the non-rotating 
case the metric corresponds to the massless multi-Taub-NUT space-time 



ds 2 



dt 



3=1 



2Nj cos 9j dipj 



dr 2 



(5.32) 



Magnetic solutions of this type (with zero dilaton field <f>) were found earlier [31]. The metric 



( 5.32 ), as well as the more general metric ( 5.26 ), may be asymptotically flat if the sum of the 



NUT charges vanishes, • Nj = 0. 



6 Type 2 solutions 



In the non-degenerate case, the matrix B, given by ( 3.21| ) where the non-collinear vectors a 
and (3 are only constrained by (3.26), depends on five parameters. One of these parameters is 
associated with rescalings of a; it is convenient to fix this scale by the condition 



detS 



(6.1) 



(the determinant of B, given by eqs. (|3.29| ) and ( [3. 32 ), is positive by virtue of eq. Q3.26D ). 
Then the four non-zero eigenvalues of B, solving the characteristic equation A| + 4 = 0, are 



±(l±i). 



(6.2) 
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The exponential e Ba may be computed using the Lagrange formula 



k=l 



Substituting (6.3) into (6.4) one obtains 



where 



e Ba = goI + 9 + B + ^{ ;]l L^ + ^D 



go = cos a cosh a, 
g± = sin a sinh a, 

g± = — (sin a cosh a ± cos a sinh a) 



(6.3) 



(6.4) 



(6.5) 



In the non-degenerate case, BPS solutions to the present system may depend only on one 
real potential, so the solution (6.4) is likely to be the most general one if deti? ^ 0. Using 
Eqs. (3.21), ( |3.25| ) and ( 3.27[ ) one obtains the exponential as a linear combination of matrices 
r a = (r a , T a i, r a i) forming a triplet of 1+2 Clifford algebras: 



where 



-Bo 



v 

V-2 



g I + v a -T a , 



-5i a A (3, 
9+& ~ 9-Ot'- 
9+P ~ 9-0- 



(6.6) 



(6.7) 



Comparing the two representations Q3.7| ) and ( |3.9| ) of M G Sp(4,R)/U(2) we obtain the 
following expressions for the potentials in terms of the components v a b of the vectors v a : 



HO 



24> _ 



r 


-1 


= 9o + + vi 2 + V2 l 


xf~ 


-1 


= vi 1 -v 2 2 , 


-vf- 


-1 


= v 1 +v 2 °, 


uf- 


-1 


= v 2 + v 1 °, 


v 2 r 


-1 


= 9o- + vi 2 - V2 l 


uvf- 


-1 


= V\ 1 +V2 2 - 



(6.8) 



Substituting Eqs. (6.6) and (6.8) one finds for the metric function / the following expression 
f^ 1 = cos a (cosh a + sinh a) + sin a (fio sinh a + fi- cosh a) , (6-9) 

where 



/i = A" 



m\ 



\4 



H± = A -1 M ± A~ 3 Re q(q 2 + 2md) 



(6.10) 



with the normalization factor A = \at A /3| 1//2 for the a, (3 of ( |4.6| ). We then obtain for the 
electric and magnetic potentials, 



v + iu = V2f (qg + + (qm + qd)A l g 1 + q(q 2 + 2md)A 2 #_j . 
The corresponding dilaton factor is 
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v 2 f 1 + cos a (cosh a + v + sinh a) — sin a (//q sinh a — v_ cosh a) 



(6.11) 
(6.12) 
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where 

v± = D ± A~ 2 Re \m{q 2 + 2md)] . (6.13) 



Points where / 1 = 0, i.e. 



cosh a -\- U-l sinh u .„ _ t . 

tan cr = : — — (6.14) 

/i- cosh a + fj,Q smh a 



correspond to spacetime singularities^]. Since f~ l enters into the spacetime metric ( |2.2| ) as 
a conformal factor of the three-space, these singularities (which are also singularities of the 
electric and magnetic potentials ( 6.11D ) are point singularities. For a one-center solution a = 



k/r, the range of r is therefore ]ri,+oo[, where n > is the largest value of r solving eq. 
( |6.14[) . The behaviour of the static (x = 0) metric near r = r% is 

ds 2 ~ c~ 2 p- 2/3 dt 2 - dp 2 - rjc 2 p 2/3 dti 2 (6.15) 

with p 2 / 3 oc (r — r±). One can show that all non-radial geodesies, as well as radial timelike 
geodesies, are deflected away from p = (r = r\). Only radial, lightlike geodesies terminate at 
the singularity p = 0, so this is a very weak singularity. Again, the corresponding string metric 
ds 2 = e 2 *ds 2 turns out to be regular. Generically the dilaton factor e 2 ^ has the same poles as 
/ so that, near r = r\, 

a 2 

ds 2 ~ t ttt dt 2 - dr 2 - rj dfl 2 . (6.16) 

(r — r\Y 

All timelike or lightlike geodesies are deflected away from the throat r = r±, except for radial 
lightlike geodesies which reach the throat after an infinite "affine time". 

The same local solution ( |6.8| ) may also be extended to the ranges r e]rfc + i, r^f, where 
Ofc = a{rk) and Ofc+i are two consecutive roots of ( |6.14 ); the corresponding spatial sections 



are compact, with two point singularities r = and r = r^+i- The more general n-center 
harmonic potential a = J2j=i kj/\ r ~ r j\ leads to an asymptotically flat/Taub-NUT solution 
with p point singularities, where p is the number of connected components of the "surface" 
a(r) = <7i, as well as compact solutions with point singularities for < a{r) < a^+i- For the 
corresponding string metrics, the point singularities are replaced by throats at infinite geodesic 
distance. 

These non-degenerate solutions may have scalar charges independent of the electromagnetic 
ones (up to the BPS condition). Let us consider some particular cases starting with type 2b 
(a, f3 non-null) which leads to surprisingly simple formulas. Here one discovers counterparts 
of the type 1 extreme Reissner-Nordstrom black holes in the dilaton-axion theory which do 
not possess dilaton charges. 

1. Electrically charged "Reissner-Nordstrom": 

/ = e 2(*-») = e ~ CT v= sina a = M/r. (6.17) 

^2cos(0-vr/4) cos(a - vr/4) ' 1 V ; 

It has Q = M, all the other charges being zero. The point singularity of the Einstein metric 
is located at r% = 4M/3tt, which is also a singularity of the electric and dilaton fields. Note 
that the non-zero dilaton field is short range. The string metric is, as discussed above, regular 
with a throat r = 4M/3ir at infinite geodesic distance. In contrast, for the type 1 extreme 
Reissner-Nordstrom solution with — D = M = Q/\^2, both the Einstein and string metrics are 
singular (the coefficient 7 in ( 5.15[ ) vanishes). 



3 In the case of the corresponding KK solutions of type 1 H], the function / _1 similarly vanishes for a 
denumerable set of values of a. These four-dimensional spacetime singularities are due to a breakdown of the 
Kaluza-Klein dimensional reduction, the five-dimensional metric being everywhere regular. 
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2. Magnetically charged "Reissner-Nordstrom" : 

/ = e -2(0+*) = ^ n= . Sin<7 7 - , a = M/r, (6.18) 

\/2cos(<t - vr/4) cos(ct - vr/4) ' v ; 

which corresponds to P = M and has a similar singularity structure. However in this case 
the product fe 2 ^ = e _2cr is non-singular, which results in the string metric ds 2 being singular 
(radial geodesies terminate at the point singularity r = 4M/3ir). 

These solutions may be subjected to 50(1, 2) x 50(2) transformations to get non-zero 
values of other charges. In particular, one can obtain massless asymptotically flat solutions 
(without NUT's). 

3. Massless electrically charged "geon" possessing a dilaton charge: 

/ = e 2 C*-*) = -= , v= , S[na j- , a = Q/r. (6.19) 

V^COS(ff-7T/4)' C0S((7 - 7T/4) ' ^' V ' 

It has D = —Q, M = N = A = P = 0. It is worth noting that the gravitational field is 
non-zero, but decays as r~ 2 . 

All the above solutions are static. 

4. Massless magnetically charged "geon" with axion charge A = P: 

f^ 1 = cos<7 cosher — sinasinhcr, \ = / (sine cosher — cos a sinh a) , 

v = u tanh a = v2/ sin a sinh cr, e 2 ^ = /(cos 2 a + sinh 2 a) (6.20) 

K = fe~ 2< ^ (sin a cosh a + cos a sinh <r) , 

where a = P/r. This solution possesses short-range electric and dilaton fields. The singularity 
corresponds to 

cotcr = tanho". (6-21) 

Note that for solutions of type 2b the only restriction on the charges is the BPS condition 
(|4.7| ). Now consider type 2a solutions (both a, (3 are null). In this case there is a second 
constraint on the charges, so that the dilaton and axion charges are functions of the other 
charges (this dependence is given by the relations (4.22), (4.23) different from those for dilaton- 
axion extreme black holes, (4.12)). 

5. Massive symmetric Reissner-Nordstrom Q = —P = without dilaton and axion 
charges: 

e _cr sino" 

\/2cos((7 - 7r/4) ' ' cos(cj - 7T/4) ' 

9 j, cosh a + e~ a sin a cos a e~ a sin a cos a — sinh a 

e 20 = 7= , k = — : = , a = M/r. (6.22) 

v2cos(cr — 7T/4) cosh cr + e °" sin a cos cr 

Note that, as in the non-dyon cases, although the dilaton and axion charges are zero, there 
still exist short-range dilaton and axion fields. 

6. Massless symmetric dyon P = —Q = Ajypl without dilaton charge: 

e°" sin a 

J = —fK ; 777' u 



y/2 cos((7 - 7r/4) ' cos(cr - 7r/4) ' 

9 j, cosh a + e CT sin a cos a sinh <r + e CT sin cr cos a .. ,„ „„ N 

e 2 ^ = 7= , «= , cr = ^/r. (6.23) 

V 2 cos (cr — 7r/4) cosh cr + e°" sin cr cos a 

This solution corresponds to an interchange a *-* (3 in the previous solution. 
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7 Discussion 



We have presented a constructive approach to the generation of BPS solutions to gravity- 
coupled field theories admitting a cr-model representation in three dimensions. This purely 
bosonic approach has the advantage of reducing all calculations to algebraic manipulations 
with matrices representing relevant cosets (no differential equations have to be solved). It 
yields solutions in a closed form readily parametrized by physical charges and automatically 
symmetric under the dualities preserving asymptotic flatness. In addition, it provides an alge- 
braic classification of solutions according to the possible algebraic types of generating matrices. 

For the particular case of dilaton-axion gravity with one vector field (which was extensively 
explored during recent years) we have listed all stationary asymptotically flat BPS solutions. 
These fall into two classes. The first class includes the well-known multi-extreme black holes 
and their rather obvious generalizations. The corresponding Einstein metrics have singular 
event horizons (null point singularities). In the case of the solutions of the new, second class, 
dilaton and axion charges are not tightly bound to electric/magnetic charges. Their Einstein 
metrics have mild timelike naked singularities. However the associated string metrics are 
generically regular for both classes of solutions. 

Both of these classes contain massless states. In the degenerate case (first class), the 
vanishing of mass can be achieved at the expense of the introduction of NUT's (although for 
multicenter configurations the total NUT charge may still be zero) . In the non-degenerate case 
(new class), there are strictly asymptoticaly flat massless solutions. 



After this work was completed, we received a preprint by Kechkin and Yurova (hep 



th/9604071 ) in which the results of [£0| and [22| are also used to construct classes of stationary 



BPS solutions to dilaton-axion gravity. However these authors incorrectly identify the coset to 
which the matrices e Ba belong as Sp(4, R) /GL(2, R). As we have seen, the correct identification 
is Sp(4, R)/(SO(l, 2) x SO(2)), the SO(2) transformations corresponding to electric-magnetic 
duality rotations. Neither the BPS solutions depending on two harmonic potentials nor our 
type 2 solutions are considered in this paper. 
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APPENDIX A: Sp(4, R) and 1+2 Clifford algebras 

The Lie algebra of Sp(4, R) consists of ten real 4x4 matrices Y satisfying the condition 

Y T J+JY = 0, (A.l) 



where J is given by ( p.13 ). For further analysis it is convenient to introduce the full set of 



sixteen independent real 4x4 matrices as tensor products of two sets of 2 x 2 matrices 

&nv = Pn ® o~v , (^,^ = 0,1,2,3), 

with 

P0 = CO = I, Pi = <Tl= Cr, pi = CT 2 = ICTy, P3 = c 3 = a z , (A.2) 
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where a x , a y , a z are the standard Pauli matrices, and it is assumed that act on 2 x 2 blocks, 
while a ^ act inside these blocks. The following multiplication rules hold 

a 2 (Xi = <7 3) 030-2 = 0-1, a 3 ai = a 2 , a\ = —a\ = a\ = l. 

In this notation, the matrices introduced in Sect. 3 read A = eo 3 , K = e 23 , while the full Lie 
algebra sp(4, R) satisfying (A.l) is 

sp(4,R) = {e 2, eio, en, e i3 , e 20 , e 2 i, e 23 , e 30 , e 3 i, e 33 } . (A.3) 

We look for a parametrization of the coset Sp(4, R)/U(2) (choosing as representatives the 
symmetric symplectic matrices M) of the type 

M = Ae Ba , A = e 3, A 2 = I, (A A) 

or, 

e Ba = AM. 

Since A G Sp(A, R)/U(2), the matrix S belongs to sp(4,i2), 

Ji?J = B T , (A.5) 
and is subject to the symmetry condition 

ABA = B T . (A.6) 

The exponential e Ba is therefore a non-symmetric matrix which can be shown to belong to the 
coset space Sp{4,R)/(SO(2) x 50(1,2)). Indeed, solving the equations ( A.6 ) and ( |A.5| ) one 
finds that B € B, with 

B = {e 2, eio, ei 3 , e 2 i, e 30 , e 33 } . (A. 7) 

The remaining generators of 5p(4, i?) form a subalgebra 

H' = {en, e 20 , e 23 , e 3 i} (A. 8) 

which can be identified with so(2) © so(l, 2), where the so(2) component is = e 23 , while the 
so(l, 2) consists of 

S a = {e 20 , -en, e 3 i} (A. 9) 

(a = 0,1,2 is an so(l, 2) index) with So = J, so that 

[E a , E b ] = 2e afe c S c , [if, E a ] = (A.10) 

(the three-dimensional Levi-Civita symbol is defined by £012 = 1, an d the indices are raised 
and lowered with the 1+2 metric rj a b = rj ab = diag(l, —1, —1)). All the elements of TL' satisfy 

AH' A = -H' T , (A.ll) 

which ensures that the similarity transformations of B generated by H' preserve the symmetry 
condition (|A.6|) . 

Inside B one can identify two sets of 1+2 Clifford algebras: 

r x a = {e 2 i, -eio, e 30 } , 

r 2 a = {e 2, e 33 , e i3 } , (A.12) 
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satisfying the anticommutation relations 

{rV r 1 b } = {r 2 a , r 2 6 } = -2 Vab i. 

Their relation to the so(2) subalgebra is clear from 



(A.13) 



k r 1 



r 1 r 2 



— 2F 
2F 1 

= 2Kr]ab, 



(A.14) 



that is, K mixes the two sets. Contrary to this, the so(l,2) does not mix them: 



s a , r b 
r 1 r 1 



2e a 6 C r 1 c 



9 c-p2 



(A.15) 



-- v\, 

Remarkably, the anticommutator of T 1 and T 2 generates a third 1+2 Clifford algebra, r°: 
{r 1 ,, T 2 b ] = -2e ah c r° c , T° a = {e 03 , e 32 , e 12 } , (A.16) 
consisting of matrices which do not belong to sp(4,R), and satisfying 

{r° a , T\} =2rj ab I. (A.17) 
This new set commutes with the so{2) £ 7i' 



K, T\ 



while for commutators with the so(l,2) E Tt' one has, similarly to (A.13), 



(A.18) 



TO tO 
1 a, 1 I 



2£ab T c, 



"2e a fe c S c . 



(A.19) 



Comparing Eqs. (A. 11), (A.13), (A.14), (A.17) one can see that the two-index 50(1,2) co- 
variant relations ( |3.24j ) hold indeed. Now, taking the commutator 



■pc -pa 



-2 V cd e ab e Z e - 2e cd e rj ab K e , 



(A.20) 



one finds that the previously introduced K is just the "time" component K = K° of the 
"vector" 

K a = {e 23 , -e i, e 22 } (A.21) 

which generates the contravariant 50(1,2) algebra dual to the "covariant" 50(1,2) algebra 
generated by the S a : 



K a , K b 

K a , T b c 
[K a , E b ] = 0. 



2e ab c K c , 



2e ab e T e c , 



(A.22) 
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The previously introduced tensor set r a & can be expressed as the tensor product of the two 
vector sets 

K a T, b = TV (A.23) 

Together with the unit matrix, the set {K a , £&, T a b} constitutes a complete basis of the space 
of real 4x4 matrices. 



APPENDIX B: Sigma model construction of the Israel- Wilson- 
Perjes solutions 



We show here how the well-known Israel-Wilson-Perjes [23] solutions of the Einstein- 
Maxwell field equations may be recovered from the general sigma model construction outlined 
in Sect. 2. 

We first introduce the familiar Ernst potentials 01 



£ = f + i X -**, $ = -j=(v + iu), (B.l) 



in terms of which the target space metric ( |3.6D reduces, for vanishing dilaton and axion fields, 

(f) = K = 0, tO 

dl 2 = —=\d£ + 2<j>d$\ 2 - -d<S>d&. (B.2) 
2/ 2 / 

This may be identified as the metric of the symmetric space SU(2,1)/S(U(2) x U(l)). A matrix 
representative of this coset may be chosen as ul 



/ 1 y/2<& i(£ -£ + 2$¥)/2\ 

M = r 1 \ -(£ + £ -2<f>$)/2 -iV2£$ . (B.3) 

\i{£ -£ -2$¥)/2 iV2£^ ££ ) 

This matrix is hermitean, M + = M, and belongs to SU(2,1), M+JM = J, detM = 1, with 

/0 —i \ 

J= | 1 . (B.4) 
\i / 

An asymptotically flat (or Taub-NUT) solution of the EM field equations depending on 
one harmonic potential a may be written in the exponential form ( |2.12| ), where the constant 
matrix 

/I ox 

A = -1 (B.5) 
\0 1/ 

commutes with J. The matrix B belongs to su(2,l), 

B + J + JB = 0, TrB = 0, (B.6) 
and is further constrained by the symmetry condition 

B + A-AB = 0. (B.7) 



These conditions also imply deti? = 0. The characteristic identity (3.32) for the matrix B then 
reduces to 

B 3 = -Tr(B 2 )B, (B.8) 
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so that the exponential representation ( 2.12j ) may be rewritten as 

M = A [1 + t~ x B sinh to + r 2 B 2 (cosh to - 1)] , (B.9) 

where t 2 = (l/2)Tv(B 2 ). 

Solutions corresponding to null geodesies in target space satisfy the further condition 

Tr(B 2 ) = 0. (B.10) 



The matrices B solving conditions ( |B.6| ), (|B.7| ) and ( B.10 ) belong to a single matrix type, with 
-B 3 = 0, B 2 0, depending (up to an overall multiplicative constant) on two real parameters 

/ I cos a —2 sin a \ 

B = -V2e^ iV2e i/3 . (B.ll) 

\— 2 sin a iy/2e~ lf} —2 cos a/ 

Two such matrices B and B' are orthogonal if 

Tt(BB') = 8[cos(a - a') - cos(/3 - /?')] = . (B.12) 

The commutator [B,B'j does not vanish (unless B' = ±B), but commutes with B and B' iff 

a ' + f3' = a + f3, (B.13) 

which also solves ( B.12 ). So the matrices B a belonging to the class ( |B.13| ) may be used to 
construct a solution (2.19) depending on several harmonic functions o a . However only two of 
these matrices are linearly independent. One can choose a basis consisting of B and of the 
matrix C defined by a 1 = a — ir/2, 0' = (3 + vr/2, e.g. 

(2 sin a —iy/2e~ l13 2 cos a \ 
-iV2e^ -y/2e ip . (B.14) 

2 cos a V2e~ if3 -2 sin a) 

This choice, such that 

BC + CB = 0, C 2 = B 2 , (B.15) 
leads to the simple form for the solution depending on two potentials 

M = Ae Ba+Cr = A[l + Bo + Ct + \b 2 {o 2 + r 2 )] . (B.16) 

Combining the two real potentials <r, r into a single complex harmonic potential 

( = o-ir, (B.17) 

we recover from ( |B.16| ) and ( [B.3D the Ernst potentials 

1 - e ia C e~ i/3 C 
£= . I , ^ , (B.18) 

1 + e ia C 1 + e ia C 

leading to the gravitational potentials 

r 1 = |l + e ia C| 2 = l + 2(o-cosa + Tsina)+<j 2 + T 2 , (B.19) 
X = 2/ (t cos a — o sin a) , (B.20) 

to be compared with the solution ( |5.4| ) for dilaton-axion gravity. Eq. ( [B.18D implies the linear 
relation between the Ernst potentials 

$ = i e -<(«+/»)(i_£) 9 (B.21) 



which is equation (24) of Israel and Wilson |23|], while their harmonic potential 2/(1 + £) — 1 
is our e m C- 
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